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In an externally applied magnetic field, ultra-pure crystals of the bilayer compound Sr3Ru207 
undergo a metamagnetic transition below a critical temperature, T*, which varies as a function 
of the angle between the magnetic field H and the Ru-0 planes. Moreover, T* approaches zero 
when H is perpendicular to the planes. This putative "metamagnetic quantum critical point", 
however, is preempted by a nematic fluid phase with order one resistive anisotropy in the ab plane. 
In a "realistic" bilayer model with moderate strength local Coulomb interactions, the existence of 
a sharp divergence of the electronic density of states near a van Hove singularity of the quasi-one- 
dimensional bands, and the presence of spin-orbit coupling results in a mean-field phase diagram 
which accounts for many of these experimentally observed phenomena. Although the spin-orbit 
coupling is not overly strong, it destroys the otherwise near perfect Fermi surface nesting and hence 
suppresses spin-density-wave (SDW) ordering. 

PACS numbers: 71.27.+a, 74.70.Pq, 71.10.Hf, 71.10.Fd, 



The Ruddlcsdcn-Popper series Sr„ + iRu„03„+i is a 
family of ruthcnate materials which exhibit a wide vari- 
ety of electronic properties ranging from unconventional 
superconductivity (n = 1) to itinerant electron ferromag- 
netism (n > 3), and have been the focus of intense re- 
search for over a decade 1 -^. Most materials in this class 
are "bad metals^" (the resistivity p(T) oc T and in excess 
of the Ioffe-Rcgcl limit) at room temperature and above; 
yet, they obey Fermi-liquid theory at low temperatures 
(typically T < 50K). The low energy electronic prop- 
erties of these materials are determined mainly by the 
electrons in the Ru ti g subspace consisting of the nearly 
degenerate d xz ,d yz ,d xy orbitals. Therefore, in addition 
to the spin and charge degrees of freedom, the orbital de- 
grees of freedom play an important role in determining 
the properties of these systems. Here, we focus on the 
bilayer (n = 2) compound Si'3Ru207 which is neither a 
superconductor nor a ferromagnet. 

Sr3Ru20ris a tetragonal material consisting of RuC>2 
planes forming bilayers which are stacked and weakly 
coupled to one another. In crystals of high purity and 
structural perfection, a metamagnetic transition^, i.e. a 
sudden and sharp rise in the magnetization with a mod- 
est increase in the applied field, is observed. While this 
transition is first-order, the transition line terminates at a 
critical point (H*,T*), where it becomes continuous (in 
analogy with a liquid-vapor transition in the pressure- 
temperature plane). However, the critical field and tem- 
perature H* and T* , depend on the angle, 9, between 
the magnetic field and the crystalline c-axis, perpendic- 
ular to the RuC>2 bilayers: H* decreases from ~ 7.8T to 
~ 5. IT , as 9 increases from to 90°, while T* drops from 
T*(9 = 90°) = 1.25K to T* ~ as 9 -> 0, that is when 



H is perpendicular to the RuC>2 planes^. Thus, it was 
proposed that this material exhibits a new type of quan- 
tum critical phenomena associated with the termination 
point of the first order line of metamagnetic transitions^. 
However, experiments involving ultra-pure, single crys- 
tal samples (with a residual resistivity less than l//f2cm) 
have shown that instead of such a "metamagnetic quan- 
tum critical point"—, there is a bifurcation of the meta- 
magnetic phase boundary, which leads to two first order 
metamagnetic transitions at closely spaced field values 
H cl ps 7.8T and H c2 ps 8.1T-2. At intermediate fields, 
H c i < H < H c2 , an electron nematic phas o 10 ' 11 appears, 
which spontaneously breaks the discrete square lattice 
rotational symmetry from C4 to C2 as inferred from the 
observation of resistive anisotropy in the ab planed 2 -. The 
nematic phase occurs in a narrow range of fields, and for 
a range of angles < 9 < 40° (the green region in Fig. 
[!}. Resistive anisotropy is also found in the blue region 
in Fig. [TJ for 55° < 9 < 90°, but it is not known whether 
this reflects the existence of a new phase. 

In the present paper, we study the microscopic origins 
of the weak metamagnetism and the accompanying ne- 
maticity. A possible microscopic route to understanding 
metamagnetism in this material was proposed by Binz 
and Sigrist 1 ^ in a model of a two dimensional band on a 
square lattice whose Fermi surface lies close to a van Hove 
(vH) singularity. Incorporating weak local Coulomb re- 
pulsion between the electrons, they showed that when 
the magnetic field tunes the Fermi surface of one spin 
species close enough to the vH singularity, there is a jump 
in magnetization. Grigera et aU^ were the first to pro- 
pose that the existence of the vH singularities might be a 
driver for nematicity accompanying metamagnetism in a 



spin-dcpcndcnt version of the Pomcranchuk distortion^ 
previously studied in the two-dimensional Hubbard and 
t-J model o 16 ' 17 . A critical insight into the problem of ne- 
maticity accompanying metamagnetism was proposed in 
a paper by H.-Y. Kee and Y-.B. Kimi&, who showed that 
additional interactions between the electrons can lead to 
new instabilities which split the metamagnetic transition, 
leading to an intermediate phase. Specifically, they con- 
sidered a model with weak quadrupolar interactions and 
showed that a mean field treatment naturally leads to a 
sequence of two transitions as a function of increasing H . 
First, at a critical field H c \, the spin-up Fermi surface re- 
connects across the vH point at one edge of the Brillouin 
zone, leading to a metamagnetic jump accompanied by 
the spontaneous breakdown of C4 symmetry (depending 
on which vH point reconnects). Then, at a higher field 
H C 2, the spin- up Fermi surface reconnects across the sec- 
ond vH point, leading to a second metamagnetic jump 
and a restoration of C4 symmetry. 

The model originally considered by Kee and Kirn^ is 
a single band model, with strongly k dependent effective 
interactions, engineered to promote nematicity. While it 
ties metamagnetism and nematicity in an ingenious way, 
it still leaves open the issue of the possible microscopic 
origins of such phenomena. In particular, it does not 
address the issue of what features of the material's elec- 
tronic structure are most important in accounting for its 
phase diagram. 

Here, we address this issue by considering a more real- 
istic model of the electronic structure of Sr 3 Ru207. Since 
both band-structure calculations and angle-resolved pho- 
tocmission (ARPES) experiments show that there are at 
least three bands at the Fermi energy per Ru atom (6 per 
bilayer), we consider a model (Section|T| with three Wan- 
nier functions, corresponding loosely to the Ru d xz , d yz , 
and d^y orbitals (Fig. [2]). The former two give rise (as 
we discuss in Section |H| to crisscrossing quasi ID bands, 
while the latter gives rise to quasi-2D bands. Each of 
these bands is, moreover, split in two by the inter-bilaycr 
hopping, which is substantial in the case of the quasi- 
1D bands. The multi-orbital band-structure also implies 
the existence of an on-site spin-orbit coupling of mod- 
erate strength. We study the ordered phases produced 
by physically reasonable local (on site) Coulomb interac- 
tions using unrestricted Hartree-Fock wavefunctions (See 
Fig. 01) 

In the multiband context, the nematic phase corre- 
sponds to a particular orbital- ordered broken symmetry 
configuration. Moreover, we will show below that transi- 
tions into such a nematic phase are naturally accompa- 
nied by metamagnetic transitions. Because the ID bands 
are closer to the vH points, and because of the stronger 
divergence of the density of states (DOS), v, at the vH 
point in ID, (y ~ \j\f~E\ we find that both nematic- 
ity and metamagnetism order are primarily driven by a 
collective reordering of these bands. In addition, a num- 
ber of other qualitative features of the experimentally 
observed phase diagram occur naturally and gcncrically 



from the mean-field solution of the present model: 

1) The interval in which the nematic phase occurs can 
be tuned to be relatively small, (H C 2 — H C {]/H C 2 "C 1. 
(See Fig. 1) 

2) There is an asymmetry to the problem, apparent 
in Fig. [5] and in the experimental data, which results 
in a monotonic decrease of the nematic order as H rises 
from H c i to H C 2, resulting in a smaller change in the 
nematicity at H C 2 than at H c \ . This feature arises in our 
model due to the underlying asymmetry in the electronic 
density of states near the vH singularity of the quasi-lD 
bands. 

3) Our model naturally accounts for why metamag- 
netism and nematic phases do not occur in the mono- 
layer ruthenate Sr2Ru04. The strong bilayer splitting 
in Sr3Ru2 07places the Fermi level in a region where the 
density of states of the ID bands has pronounced posi- 
tive curvature, satisfying the requirement of the Landau 
theory of a weakly first order metamagnetic transition 19 . 

4) Even a moderate spin orbit coupling, A s . . ~ 0.2<, 
consistent with band structure estimate a 20 ' 21 , produces 
an order 1 decrease of the critical fields as 9 varies from 
to 90°. (See Fig. E) 

5) The near-perfect nesting of the quasi-lD bands ac- 
counts for the most prominent peaks in the low energy 
magnetic structure factor observed in neutron scattering 
experiments at H = 22 i 23 . However, the spin-orbit cou- 
pling is remarkably efficient at spoiling this nesting, thus 
plausibly explaining why spin density wave (SDW) order 
does not actually materialize (at least for H = 0) 22 i 23 . 

Clearly, there are many aspects of the physics that are 
more subtle, and cannot be addressed, even qualitatively, 
at mean field level. There is, after all, considerable evi- 
dence of the effects of strong quantum fluctuations asso- 
ciated with the narrowly preempted metamagnetic quan- 
tum critical point. We will return to these shortcomings 
at the end of the paper. 

I. THE MODEL 

We have studied a simple tight-binding model of a sin- 
gle RUO4 bilayer, with terms organized according to a 
hierarchy of scales: 

H = H l +H s . .+H 2 + ... (1) 

where H\ contains the largest terms, which involve the 
most direct 7r-overlaps between Ru d-orbitals on nearest- 
neighbor sites: the largest hoppings between two neigh- 
boring identical d aa > orbitals (a, a' = x,y, z, a ^ a') are 
along the crystalline a and a' directions. These hoppings 
in turn make use of the intervening oxygen p-orbitals. 
The single bilayer approximation is a good one because 
transport measurements confirm the existence of highly 
two-dimensional transport, and therefore weak bilayer- 
bilayer coupling. Hi also includes the onsite Coulomb 
repulsion terms between two electrons on the same or- 
bital (U), as well as between two electrons in different 
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FIG. 1: (Color online) The experimentally determined low 
temperature phase diagram of SrsRu2 07in the field - angle 
plane, based on resistivity and magnetic susceptibility mea- 
surements on shape-unbiased octagonal crystals at 100 mK. 
For the resistivity measurements, both the in-plane field com- 
ponent and the current are along either the crystallographic 
a or b axes. The shaded regions are those in which a resistive 
anisotropy is observed. The solid black lines represent first- 
order phase transitions as determined by a sharp dissipative 
peak in a measurement of the imaginary component of the a.c. 
susceptibility [see discussion in Grigera et aZ.— ]. We observed 
no dissipative peak at the boundaries of the blue region, but 
cannot rule out its existence beyond our resolution. 




FIG. 2: (Color online) The Ru d xy ,d yz ,d x ~ orbitals in a sin- 
gle layer of Sr3Ru207in the ab plane, rr-overlaps between 
two like orbitals along nearest-neighbor bonds are mediated 
by the intervening oxygen p-orbitals. For two identical dij 
orbitals (i ^ j, and i, j — x, y, z), the hopping along the i di- 
rection is strong and is mediated by the oxygen pj orbital and 
vice-versa. 5-overlaps between two identical d-orbitals along 
a nearest neighbor bond do not make use of the oxygens, 
and are therefore much weaker in comparison. Furthermore, 
all nearest-neighbor hopping between two distinct d-orbitals 
vanish by symmetry. 



orbitals (V). H s o , captures the effect of onsite spin-orbit 
coupling. H2 represents the kinetic energy terms due to 
weaker 5-overlaps between the orbitals (e.g. along the z 
direction for the d xy orbitals). Still smaller terms, some 
of which we will mention below, are represented by the 
ellipsis. 

Let d} ^ ^ create an electron with spin polarization 

a at horizontal position i?, in layer A = ±1, and in or- 
bital a = x,y,z corresponding to the d xz , d yz , or d^ 
orbital respectively. In order to emphasize the undcrly- 



FIG. 3: (a) A "realistic" tight-binding Fermi surface of 
Sr3Ru207taking into account the Kinetic energy terms in Hi, 
Hi as well as H a . .. Here, t' = O.lt and A s . . = 0.2t. (b) 
Fermi-surface derived from an idealized model taking into ac- 
count only the single particle terms in Hi. In the nematic 
phase, one of the "spin-up" bonding bands crosses the vH 
point at either (it, 0), or (0, 7r) depending on the sign of the 
nematic order parameter (to be defined below). Quotes are 
placed around "spin-up" because in the presence of spin-orbit 
coupling, the magnetic field lifts the degeneracy of Kramers' 
doublets, and in this case, it is the appropriate pseudo-spin 
band which crosses the vH point. 

ing symmetries of the Hamiltonian, we define a spinor 
field 

*L(£>A) = < iCTji?)A , a = x,y,z. (2) 
In terms of these, and n X aM = £ CT d [^^ x d a ,a.,R.x-- 
#i = -t^2 [* f (i?,A) (T x + T z ^j $(R + x,\) 

A, ft 

+^(R, A) (fy + f z ^j %{R + y, A) 

+ -&(R, -A) (f x + T s ) $(R, A) + H.C. 

H — 7 7i, ^ H > 7i, gn, , s 

A, a, ft. A,a^a'ft 

-h-J2^(R,^)(l + s)^(r,x) (3) 

ft,A 

The contraction over the spinor subscript indices is im- 
plied and we have defined the following matrices: 
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where £ are the orbital angular momenta projected onto 
the t2 g states (an explicit form of these are given in Sec- 
tion [TTT] ), and t are the Pauli matrices. The final term 
above is the Zeeman coupling to an external field. We 
note that in the tig subspace, the angular momentum is 
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only partially quenched, since it is possible, for example, 
to form linear combinations of orbitals d xz ± id yz which 
are eigenstates of £ z . Therefore, the external magnetic 
field will couple both to L and S. 

The on-site spin-orbit coupling Hamiltonian is 

H s .o. = \ s .o. A) (L ■ S) A) (5) 

ft, A 

and the smaller nearest-neighbor couplings arc 



H 2 = -t'J2^(R,X)f x ^(R + y,X) 

R,X 

+f t (i?, \)f y <f{R + x,\) 

+ -&(R, \)T z <t>{R, -A) + H.C. 



(6) 



So far, in writing the above terms, we have assumed lo- 
cally perfect octahedral symmetry with the result, for ex- 
ample, that the hopping matrix element, t, between d xz 
orbitals on sites separated by one lattice constant in the 
x-direction, or on equivalent sites in neighboring planes 
of a bilayer are equal to each other. The actual material 
deviates slightly from this ideal symmetry^; small terms 
that break this symmetry, as well as further-range hop- 
ping terms and interactions, are all represented schemat- 
ically by the ellipsis in Eq. [1] and will not be considered 
explicitly here. 

The Fermi surface of our model taking into account 
the kinetic terms of H\ , H 2 and also the spin-orbit term 
H s o . is shown in Fig. [3^. It is useful to contemplate its 
relation to the simpler Fermi surface (shown in Fig. [3j} ) 
obtained by setting the couplings in H2, and H s o , to 0. 
In this simpler model, hoppings only along the strongest 
bonds of each orbital are taken into account, and conse- 
quently, there arc perfect one-dimensional bands which 
are split by the bilayer hopping and form the straight 
patches of the Fermi surface. These Fermi sheets are 
purely d XZl d yz in orbital character. Moreover, the 2D 
bands which come from the d xy orbital are degenerate 
at this level of approximation. Were we to include the 
effects of the smaller hopping terms in H2 , a small split- 
ting of the 2D bands, and a slight warping of the ID 
bands would result, but there still would not be any mix- 
ing between bands, and hence there would be multiple 
points at which two pieces of the Fermi surface would 
cross one another. When spin-orbit coupling is included, 
these degeneracies are lifted, as is evident from Fig. (3Ju 
Therefore, i? s . . has the important qualitative effect of 
changing the Fermi surface topology, even if its magni- 
tude is small. We will study the consequences of includ- 
ing H s o in Section HTT1 

The d xy bands play little active role in the physics dis- 
cussed below, so we will for the most part ignore these 
bands altogether, treating them instead simply as a par- 
ticle reservoir which allows us to perform calculations 
involving the remaining bands at constant chemical po- 
tential rather than at constant density. 
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FIG. 4: (Color online) Mean-field T — phase diagram in the 
U — V plane in the absence of spin-orbit coupling for fi = 0.81, 
at zero applied field (a), and in finite field (b), taken here to be 
\H\ — 0.04t. In zero field (a), there are four distinct phases: 
a paramagnetic region (I), a nematic region (shown in blue) 
(II), a region in which the ferromagnet and nematic spin ne- 
matic phases are degenerate, but do not coexist (III), and a 
region in which the nematic, ferromagnet, and nematic spin 
nematic phases coexist (purple region) (IV) . All phase bound- 
aries correspond to first order transitions. In a finite field (b), 
all regions acquire a non-zero magnetic moment. Regions I 
and III do not have nematic order, whereas region II contains 
nematic order. In this case, phase boundaries correspond to 
metamagnetic transitions. For U ~ V chosen such that at 
zero field, the parameters lie in region I in (a) above, a field 
sweep moves the phase boundaries closer to the origin so that 
the system traverses regions I — > II — » III in (b) as the field 
increases. In this way, the system exhibits a nematic phase 
precisely between two metamagnetic transitions. 



II. MEAN-FIELD THEORY 

We first consider an idealized model of a bilayer con- 
sisting of 2 orbitals (d xz , d yz ) in each layer. For the mo- 
ment, we neglect the curvature of the Fermi surface and 
any mixing between the two orbitals. In this limit, the 
orbitals form perfect one-dimensional bands that are bi- 
layer split due to the hopping along the c-axis. In the 
presence of these multiple orbital degrees of freedom, it 
is natural to consider spin and orbital ordered broken 
symmetry states, and so we define the following sets of 
collective variables: 

N\ = ^(*t(i?,A) (f x + TA *(i?,A)) 



M X = J2{&(R,\) 



S [T x + T v 



A)) 



N° x 



Y^(^{R, A) (f x - T y ) *(i?, A)) 



]>>t(i?,A) 



S(T X - T v 



*(ii,A)) (7) 
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FIG. 5: (Color online) Phase diagram at zero temperature 
in the U — H plane with U — V in the absence of spin-orbit 
coupling for /i = 0.81. For each U — V, as the field is in- 
creased, two metamagnetic transitions occur and are denoted 
Hex (dashed line) and H C 2 (solid line). For H c i < H < H C 2 a 
nematic phase (blue region) occurs, corresponds to region II 
in Fig. 2p, and collapses onto region IV in Fig. in the zero 
field limit. Although we have chosen U — V in this figure, 
we emphasize that a phase diagram with the same topology 
occurs if we set U = aV, for a range of a, so long as the line 
U = aV crosses each of the three regions shown in Fig. [4b. 
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FIG. 6: (Color online) (a) Magnetization and (b) nematicity 
as a function of applied field for U = V = 2.5t, fj, = 0.81t 
in the absence of spin-orbit coupling. The nematic phase is 
found to lie precisely between the two metamagnetic transi- 
tions. 

These represent, respectively, the total electron density 
of the one-dimensional bands, N, their overall magneti- 
zation, M , the nematic order parameter, N°, that rep- 
resents the difference in charge density in the two bands, 
and lastly, we also define N s , which represents the differ- 
ence in moment in each band. We call this latter quan- 
tity the order parameter for the "nematic spin nematic" 
phased. Whereas M breaks time-reversal (T) and SU (2) 
symmetry, and N° breaks the lattice C4 rotation symme- 
try (R), the nematic-spin-nematic order breaks SU(2), 
(T) and (R); however, the product (TR) is not broken in 
this phase. In the presence of an externally applied mag- 
netic field, there is no distinction between the nematic 
and the nematic-spin-nematic phases. 



The mean-field phase diagram of our model is pre- 
sented in Figure 2] The details of the calculations are 
presented in the appendix. The zero-field phase diagram 

(a) consists of 4 phases: I: a paramagnetic phase, II: 
a nematic phase, III: a phase in which the ferromagnet 
and nematic-spin-nematic are degenerate but do not co- 
exist, and IV: a coexistence phase of nematic, nematic- 
spin-nematic, and ferromagnetism. All phase-boundaries 
shown here are first-order transitions. When the mag- 
netic field is non-zero (b), the ferromagnetic phase has 
a lower free energy than the nematic-spin-nematic phase 
and this degeneracy is lifted. 

In a non-zero field, all phase boundaries become meta- 
magnetic transitions. Two qualitatively different features 
arise in this case. First, there are only two types of 
phases: either magnetic order is present alone (I, III), 
or magnetic order coexists with nematic order (II). In 

(b) , the only difference between regions I and III is the 
size of the moment (larger in region III). Note in particu- 
lar, that with increasing magnetic field, phase III moves 
in to smaller values of U and V. Thus, without requir- 
ing any fine-tuning, it is possible for the system studied 
here to undergo a sequence of metamagnetic transitions, 
starting with a paramagnetic phase, ending with a fer- 
romagnetic phase, with nematic order sandwiched in be- 
tween, as shown in Fig. [5] The requirements for such a 
phenomenon to occur in our model is that U ~ V. For 
U 3> V, for instance, the effect of applying a magnetic 
field is to induce ferromagnetism alone. By contrast, 
when U ~ V, this ferromagnetic phase is preempted at 
low fields by a nematic phase. Figure [HI shows the mag- 
netization and nematic order as a function of field for a 
particular value of U = V = 2.5t. Here it is clearly seen 
that as the field sweeps across the first metamagnetic 
transition, nematic order develops; when the field is in- 
creased further, a second metamagnetic transition occurs 
which destroys the nematic phase. 

In Fig. \7\ the single quasi-particle density of states is 
plotted in the nematic phase. There is a double peak near 
the Fermi level, whose splitting corresponds to the size 
of the nematic gap. Also shown in Fig. [7] is the density 
of states at a slightly lower field, where the system is just 
about to enter the nematic phase. We note that there is 
a reduction of the quasiparticle density of states at the 
Fermi level in the nematic phase; this in turn implies that 
the entropy of the nematic phase at low temperatures 
within mean-field theory is lower than the neighboring 
isotropic phases. This is also seen in Fig. [5] from the 
curvature of the phase boundaries of the nematic phase 
at finite temperatures. We highlight here that this result 
from mean-field theory disagrees with the experimental 
observations reported in Ref J^, where the nematic phase 
boundaries "fan" outward, implying that in actuality, the 
nematic phase has a higher entropy than the neighboring 
isotropic phases. We believe that this effect must stem 
from fluctuations (both thermal and quantum) due to the 
presence of nematic domains. 
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FIG. 7: Single quasiparticle density of states (arbitrary units) 
in the nematic phase, H = 0.04i, U = V = 2.5t, fx = 0.811 
The Fermi energy is at E = 0. The dashed line corresponds 
to the quasiparticle density of states at H = 0.02i, there 
the nematic phase does not arise. Since the nematic phase 
has a lower density of states, it is expected that in mean-field 
theory, it will also have a lower entropy at temperatures small 
compared to \i. 
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FIG. 8: Phase diagram of the system without spin-orbit cou- 
pling showing the finite temperature boundaries of the ne- 
matic phase. First-order transitions (solid curve) at low tem- 
peratures give way to continuous transitions (dashed curve) at 
higher temperatures as the vH singularities get smoothed out. 
Finite temperature metamagnetic cross-overs away from the 
nematic phase are not shown. Here, U — V — 2.5t, /j, = 0.81t 
and A = 0.2*. 



III. EFFECT OF SPIN-ORBIT COUPLING: 
ANGLE-DEPENDENT METANEMATIC 
TRANSITIONS 

The results above were based on the assumption that 
spin SU (2) symmetry is preserved in the system. How- 
ever, the experimental phase diagram^ of Sr3Ru 2 C>7(Fig. 
[I]) exhibits an 0(1) anisotropy in the critical field at 
which a metamagnetism occurs: the critical field with 
H along the c axis is about 1.6 times larger than the case 
when it points in the ab plane. A natural way to account 



for this anisotropy is to consider the effect of spin-orbit 
coupling in the material. Since the material has an inver- 
sion center, k-dependent spin-orbit interactions, such as 
the Rashba coupling, are forbidden, and we proceed by 
considering the effect of an atomic spin-orbit coupling, 
given in equation [5J While the precise magnitude of 
the spin-orbit coupling constant X s . , for this material 
is unclear, LDA calculations of the monoloycr ruthenate 
compound S^RuO,}, suggest that A S . D . is approximately 
10 percent of E f 20 ' 21 . Indeed, in recent angle-resolved 
photoemission (ARPES) studies of Sr 3 Ru207, LDA cal- 
culations had to employ spin-orbit coupling in order to 
fit properly to the Fermi surface^. Since this spin-orbit 
coupling term is onsitc, it seems reasonable to expect 
that similar values also hold for the bilayer compound. 

When such spin-orbit coupling terms are included in 
the Hamiltonian, the d xy band must necessarily be taken 
into account. However, since the spin-orbit coupling is 
smaller than the crystal field splitting, we can treat it 
as a perturbation, and project the angular momentum 
operator L onto the *2g subspace. To be more explicit, a 
valid choice of the orbital angular momentum operators 
projected onto the *2g manifold is 









£ z = -i[£ x ,£ y ] 



(8) 



In terms of these operators, the spin-orbit coupling in 
each layer of our system is represented by the 6x6 matrix: 



L S 



i x - i£y 



+ 



(9) 



Having in mind the framework of a minimal model 
which captures the essential features of the experimental 
phase diagram, we first treat the quasi-2D band as a free 
electron system and neglect the effect of Coulomb inter- 
actions. We have checked that the inclusion of a Hubbard 
repulsion on this band produces no qualitative changes 
to the results reported here. Furthermore, we neglect 
the terms in H2 as before. This way, the mean-field or- 
der parameters in the presence of spin-orbit coupling arc 
identical to those in Eq. [Jj and the spin quantization 
axis is defined to be in the direction of the applied field. 
However, the spin-orbit coupling term is sensitive to the 
orientation of the magnetic field relative to the c-axis 
since the orbital angular momentum operators above are 
defined with respect to the crystalline axis of the system. 
More explicitly, in the tilted field, the spin-orbit term is 
modified as 



H so {6) = \ S . L-S(6) 
.3 ■ n 



S(6) 



exp 



9)Sexp[i^d), (10) 



where 9 is the angle of the applied field relative to the 
c-axis and n is a unit vector either in the crystalline a or 
b direction. 
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FIG. 9: (Color online) Angle-dependent metanematic phase 
boundaries in the presence of spin-orbit coupling showing an 
O(l) anisotropy in the critical fields. When the field is tilted 
towards the ab plane, the C4 point group symmetry is au- 
tomatically broken. In region I at low fields, the moment is 
small and there is also a small nematic phase present for all 
6 > 0. The first metamagnetic transition into region II gives 
rise to an accompanying jump in nematicity in addition to the 
magnetization. Finally, the second metamagnetic transition 
into region III results in a discontinuous decrease in the size 
of the nematic order parameter as well as a jump in magne- 
tization. Here, U = V = 2.5i, fi = 0.81£ and X s . . = 0.2t. 



In Fig. [5J we show the phase diagram in the H — 6 
plane, keeping U = V fixed and A s . . = 0.2t. As the field 
is tilted towards the ab plane, we see that the critical field 
at which the metamagnetic transition occurs is smaller in 
magnitude. For all angles, there are two metamagnetic 
transitions with a nematic phase in between. When the 
field is tilted away from the c-axis, the crystalline C4 
symmetry is broken and we distinguish a phase in which 
the nematic order jumps from a small value (region I) 
to a large value (region II) and jumps back to a small 
value (region III) . Thus the tilted field shows the nematic 
analog of metamagnetic transitions, and we refer to these 
transitions here as "metanematic" transitions 2 ^. We note 
that although we find metamagnetic transitions with the 
same general angle-dependence as seen in experiments, 
the nematic phase in our model does not have the cor- 
rect topology in the field-angle plane (compare Fig. [9]and 
Fig. [1} . Instead, we find that if a nematic phase occurs, it 
remains present for all field orientations. Given the struc- 
ture of the phases shown in Fig. 2}d, it does not appear 
impossible that the correct topology could emerge under 
appropriate circumstances, but in the present model, this 
would require angle-dependent interactions. While such 
terms are unphysical at the bare microscopic level, they 
could arise at the effective level due to fluctuation effects 
that have not been considered here. 

The ordering tendencies which we have considered so 
far do not break lattice translational symmetry. How- 
ever, inelastic neutron scattering experiment s 22 ! 23 have 
provided evidence that there are substantial incommen- 



surate spin fluctuations in this system, although static 
spin-density-wave (SDW) order has not been observed 
at zero external field. The incommensurate spin fluctu- 
ations in this system occur primarily due to the partial 
nesting of the Fermi surface. Therefore, we must check 
whether such finite q ordering tendencies in this material 
are favored over the uniform nematic order proposed in 
this paper. To do this, we have computed the generalized 
one-loop spin susceptibility 

pp' a(3^8 

( T rdi pa {T)d tp+ql3 (t)4 p ' 7 (0)d bp >- qS (0)) 

(11) 

in the presence of spin-orbit coupling and magnetic 
field. The electron Hamiltonian is a 12 x 12 matrix, 
and the electron propagators are spin and orbital depen- 
dent which in turn makes the above spin susceptibility a 
36 x 36 matrix. Using the random phase approximation 
(RPA), we find that an instability towards the formation 
of SDW order occurs in the system (taking U — V) when 
U exceeds a critical strength U c . s d w which satisfies 

^Max(cig[ X ]) = l (12) 

and for the band-structure parameters we have been us- 
ing (A s . . = 0.2i, t' = 0, fj, = 0.81t) the maximal eigen- 
value of the susceptibility matrix is obtained at the wave- 
vector q = (0.27, 0.27) it where the spin susceptibility 
X zz (q) obtains its largest value. From this mean- field es- 
timation, we find that U c ^ s d w ~ 3.1i, which is also close 
to the critical coupling required for Fcrromagnctism, as 
seen in Figure [4j 

Thus, in addition to accounting for the angle- 
dependent metamagnetic transitions observed in this sys- 
tem, the inclusion of spin-orbit coupling also naturally 
explains why static incommensurate SDW order does not 
occur despite the presence of incommensurate spin fluc- 
tuations in this system. 



IV. DISCUSSION 

We have shown in this paper that the remarkable low 
temperature properties of SraR^Oycan be understood 
as being a consequence of orbital ordering of the quasi- 
1D bands. We have been able to account for most of 
the gross features of the experiments from a simple mi- 
croscopic model. Our results are based on a mean-field 
solution of a system with strong electron interactions and 
a justification for focusing solely on this approach has not 
been provided. Although metamagnetism and nematic- 
ity generally arise as strong coupling effects in metals and 
are therefore difficult to treat in a controlled theoretical 
fashion, the divergence of the density of states v at the 
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vH point insures that the Stoncr criterion can be satis- 
fied even for weak interactions, which are much less than 
the bandwidth. Thus, by the notion of adiabatic conti- 
nuity, it may be legitimate to treat this problem from a 
weak-coupling standpoint. 

Previous attempts to explain metamagnetism and ne- 
maticity in this syste m 18 ' 28 ' 29 were based on the assump- 
tion that the quasi-2D bands drive the nematic and meta- 
magnetic transitions. By contrast, we point out here that 
it is much more natural to think instead that the quasi- 
1D bands are responsible for the transitions. We note, 
from a symmetry standpoint, that the quasi-lD bands 
form a two-fold representation of the C4 rotation symme- 
try of this system: thus, the nematic phase which breaks 
C4 rotation symmetry corresponds to an orbital order- 
ing among these bands. Furthermore, it is known from 
experiment that the monolayer ruthenate Si'2Ru04 does 
not exhibit metamagnetism for magnetic fields upto 30 
Tesla. The primary difference in the electronic structure 
of the monolayer vs. the bilayer compound is the large 
bilayer splitting in the latter. We have shown here that 
while the quasi-2D bands are only weakly affected by the 
bilayer splitting, the quasi-lD bands are rather strongly 
affected by it. Thus, the experimental differences be- 
tween the monolayer and bilayer ruthenate compounds 
as well as symmetry considerations lead us to propose 
that quasi-lD bands in this system are primarily respon- 
sible for the rich phase diagram of the bilayer ruthenate. 
Although we propose a different microscopic origin for 
the Fermi surface distortion from that used in previous 
work on the subject, we note that our picture is still one 
of a weak-coupling Pomcranchuk type. As such it auto- 
matically retains the attractive feature, common to any 
weak anisotropic distortion of electronic structure in the 
vicinity of Ep, that it naturally predicts a sensitivity to 
impurity scattering^. This in turn matches one of the 
key experimentally determined characteristics of the be- 
havior that we set out to explain. 

Along with the focus on the ID bands comes the ex- 
istence of a nesting vector, 2kp, and the issue of SDW 
order. At mean-field level, the preferred SDW instabil- 
ity for a pair of orthogonal quasi- ID bands leads^I ei- 
ther to bidirectional order with ordering vectors (2kp, 7r) 
and (71", 2kp), or unidirectional order with ordering vector 
(2fcp, 2kp). However, inelastic neutron scattering studies 
have found that the most dominant peaks are observed 
in the (2kp,0) direction^ 2 .. Nevertheless, we have shown 
that even moderate spin-orbit coupling reduces the nest- 
ing to the extent that the mean field tendency to nematic 
and SDW ordering are comparably strong. Moreover, we 
expect fluctuation effects to reduce the ordering tendency 
of an incommensurate SDW (which has a gapless sliding 
mode and two nearly gapless spin-wave modes) relative 
to the Ising-likc nematic order, especially given the quasi- 
2D structure of this material, which makes fluctuational 
corrections all the more significant. 

However, there are still two qualitatively important 
features of the experimental data that are not well ac- 



counted for by any mean- field treatment we know of: 1) 
The intermediate nematic phase has been shown, in ex- 
periment, to have higher entropy at low temperatures 
than either of the adjacent disordered phases. In con- 
trast, at mean-field level, the ordered phase always has a 
lower Fermi surface density of states than the proximate 
disordered phases, and so has smaller low temperature 
entropy. 2) Doping studies suggest that at least part of 
the peaks in the density of states that are thought to drive 
the metamagnetism are somehow locked to the Fermi en- 
ergy, rather than being purely features of an underlying 
rigid band-structure. These features we believe are sig- 
natures of strong-coupling effects, and can be accounted 
for by a more sophisticated theory in which local nematic 
order is present over a broad range of magnetic fields (in- 
cluding B = 0), but only propagates to long distances in 
the narrow range of B in which macroscopic anisotropics 
are observed. 

More generally, it is clear that fluctuation effects play 
a significant role in the physics. Even though the puta- 
tive metamagnetic quantum critical point is preempted 
by the nematic phase, the observed fluctuational phe- 
nomena that lead to the conjectured critical point in the 
first place remain real and dramatic. It is likely that 
they reflect the existence of a "nearby" quantum critical 
point, even if it is not actually observed. In this context, 
mean-field results, of the sort discussed in the present pa- 
per, should be adopted with caution. At the very least, 
the effective parameters that enter our model must be 
reinterpreted as strongly renormalizcd effective parame- 
ters, given that the observed bandwidthsH are narrower 
by a factor of order 10-100 than the bandwidths found 
in LDA calculations^. We defer the fascinating study of 
fluctuation effects to a future publication. 

Finally, we make a comment about the observed re- 
sistivity anisotropy in this material when H || c. While 
we have presented a symmetry argument for why such a 
resistivity anisotropy ought to be present, we have not 
explicitly computed the resistivity anisotropy from our 
model. Indeed, any small nematic Fermi surface distor- 
tion like that discussed here is unlikely to produce a large 
effect on the low-temperature resistivity due to the dis- 
tortion alone. The reason for this, is that the nematic 
phase arises from a discrepancy in electrons which are 
close to a vH singularity; these, in turn, have a very 
small characteristic velocity and hence would contribute 
most weakly to transport signatures. The distortion of 
the Fermi surface is more likely to change the transport 
properties by being a source of domains and domain-wall 
scattering, something which is beyond the scope of the 
calculations that we report. We note further that in the 
experiments of RefJ^, when the field is tilted towards the 
ab plane, both the average resistivity and the anisotropy 
rapidly decrease with angle. In our picture, we imagine 
that as the field is tilted, nematic domains get aligned 
and scattering is therefore considerably reduced. Thus, 
ironically, precisely the same signal which was used to 
detect the nematic fluid would prove to be useless deep 
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in the nematic phase, when the system forms a single 
macroscopic nematic order. A more quantitative theory 
of the physics discussed here will be presented elsewhere. 
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APPENDIX: METAMAGNETISM AND 
NEMATICITY OF QUASI-1D BANDS 

Here, we present, for the sake of clarity, the derivation 
of the Mean-field equations which were used to deduce 
the phase diagram in Fig. [4] The mean-field order pa- 
rameters are defined via 



-(N x 
4 V 

1 



in 



y,X.R,al 



= - (N x - N x 
4 V o 



aM x + aN x ) 
aM x -aN x ) (A.l) 



After decomposing the interactions in terms of the above 
expectation values, we arrive at the following one-particle 
Hamiltonian: 



H(k) 



H x (k) 

Hy(k) 



where 



H a {k) = -2icosfc a i 4x4 + 
( 5/j$ - h -t 
h 







(-1) 











and a = x,y. We have also defined the quantities 



V - (N x - o~M x - oN x ) - ^N x 

j {-N x - oM x + aN x ) + | N x (A.3) 



Thus, by neglecting the hybridization between the x, y 
orbitals, the Mean- field Hamiltonian takes the Block- 
diagonal form above and the quasiparticle bonding and 
anti-bonding energies are easily obtained: 



—2tcosk x — 
± 



x ah + 1 + 



^ 2 + i(^-^; (J 1) )' 



e y ±a = -21 cos k y - ah+ ^ Um^I + fyy, 



(A.4) 



Due to the ID band dispersion, it is possible to obtain 
the density of states and the grand-canonical free energy 
density analytically by summing over the 8 quasi-particlc 
bands of this model. At zero temperature, 



F = F (M X 1 N X 1 N X )+ f {E-n)u[E)dE 

J —OO 



(A.2) F Q (M X , N X ,N X ) = ^ [(^) 2 + ( M> 



V 
1 



N 



Nt 
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The phase diagram is then obtained by minimizing the 
above free energy with respect to the order parameters. 
We note that even upon including curvature effects in 
the band-structure, the grand canonical free energy den- 
sity can be expressed analytically in terms of complete 



Elliptic Integrals. However, for the sake of simplicity, we 
don't include such terms here. While such Fermi terms 
do modify the precise location of phase boundaries, they 
have no qualitative effect on the physics. 
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